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ABSTRACT 

The Griffith-Woodley Technique (GWT) is an approach to estimating precipitation using 
infrared observations of clouds from geosynchronous satellites. It is examined in three ways: an 
analysis of the terms in the GWT equations; a case study of infrared imagery portraying convec- 
tive development over Florida; and the comparison of a simplified equation set and resultant rain 
map to results using the GWT. The objective is to determine the dominant factors in the calcula- 
tion of GWT rain estimates. 

Analysis of a single day’s convection over Florida produced a number of significant insights 
into various terms in the GWT rainfall equations. Due to the definition of clouds by a threshold 
isotherm the majority of clouds on this day did not go through an idealized life cycle before 
losing their identity through merger, splitting, etc. As a result, 85% of the clouds had a defined 
life of 0.5 or 1 h. For these clouds the terms in the GWT which are dependent on cloud life 
history become essentially constant. The empirically derived ratio of radar echo area to cloud' 
area is given a singular value (0.02) for 43% of the sample, while the rainrate term is 20.7 mm IT* 
for 61% of the sample. For 55% of the sampled clouds the temperature weighting term is identi- 
cally 1 .0. Cloud area itself is highly correlated (r = 0.88) with GWT-computed rain volume. An 
important, discriminating parameter in the GWT is the temperature defining the coldest 10% 
cloud area. 

The analysis further shows that the two dominant parameters in rainfall estimation are the 
existence of cold cloud and the duration of cloud over a point. This leads us to conclude that the 
GWT is unnecessarily complicated for use in estimating daily rainfall over iarge areas. Simplifying 
assumptions are made to the GWT such that the resultant equations are independent of cloud life 
history, cloud area, and grid square area. Application of a simple algorithm incorporating these 
assumptions leads to a daily rainfall pattern very similar to that calculated from the GWT for one 
case study day. We present an in-depth analysis of the storm which caused the only significant, 
difference between the two isohyetal maps. 



RAiN ESTIMATION FROM SATELLITES: AN 
EXAMINATION OF THE GR1FFITH-WOODLEY TECHNIQUE 

1. INTRODUCTION 

The recent Workshop on Precipitation Measurements from Space (1981 ) at NASA/Goddard 
Space Flight Center stressed the need for spaceborne measurements of precipitation across the 
spectrum of scales from Hood warnings to global climate systems. The workshop participants dis- 
cussed a variety of visible and infrared (IR) methods for observing cloud properties and using them 
as proxy variables for precipitation. One such method is the Griffith- Woodley Technique (GWT), 
an approach to estimating precipitation from observations of clouds in the 1 lfj window channel 
from geosynchronous satellites (Griffith et al„ 1976). Briefly, the GWT: 

1) defines and tracks clouds by the outline of the 253K isotherm; 

2) estimates a cloud rain volume (R ) through an algorithm that: 

a) relates the rain area (echo) to the life history of the cloud area; 

b) relates the rain rate to the life history of the rain area (echo); 

c) enhances the R v by a weighting term based on the distribution of cloud-top temperature; 

3) produces a daily rain map by apportioning the calculated R v to grid squares by concentrating 
the rain in the coldest portions of the cloud. 

The purpose of this paper is to examine the GWT in three ways: an analysis of the magnitude 
and variability of the terms in the GWT equations; an examination of infrared (IR) imagery in terms 
of GWT-defined clouds and their life cycles; and thirdly, the application of a simplified equation 
set and the resultant comparison to the GWT results. The objective is not to improve the GWT, nor 
to propose a new technique. Rather, wc wish to determine what are the dominant factors in the 
production of GWT rain estimates. 
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At the core of the original GWT are empirical relationships between the temporal history of 
satellite-determined visible cloud area and gage/radar measured rainfall, over southern Florida 
(Griffith etal., 1976). These relationships were redefined as functions of cloud-top temperature 
(CTT) when digital IR data became available. Rainfall estimates were subsequently determined 
for southern Florida, tropical regions in South America, and for selected Atlantic hurricanes (Grif- 
fith et al., 1 978). A scheme to apportion the total rain volume of cloud segments were constructed 
(Woodley et al., 1980) and modified (Augustine et al., I98i ) and isohyetal maps of GATE B-scale 
rainfall, were generated. The GWT has also been applied to the U.S. High Plains (Griffith et al., 
1981). These results were adjusted in one of two ways: by gage/satellite comparisons for a small 
area; or by noting differences in the local sounding compared to a mean tropical sounding. This 
latter adjustment factor was based on output from the Simpson-Wiggert (1969) one-dimensional 
cumulus cloud model. Recently, Meitin et al., (1981) have used the GWT to construct isohyetal 
maps for 51 days during the Florida Area Cumulus Experiment (FACE-2). 

That rainfall and cloud parameters are highly correlated, especially over large time and space 
scales, has been shown by several investigators. Arkin (1979), for the GATE B-array (10 s km 2 ), 
showed correlation coefficients as high as 0.88 between 6 h accumulated rainfall and the fraction of 
the array covered by cloud higher than 10 km (based on an IR threshold CTT). This compares to 
the coefficient of 0.87 found by Woodley et al. (1980) between rain estimated by the time and 
temperature dependent GWT and verification rain from radar. Garcia (1981) estimated GATE 
rainfall using the simple Kilonsky-Ramage (1976) regression equation in conjunction with one 
satellite image per Jay. Cumulative GATE phase estimates (17-20 day period) were highly corre- 
lated with GWT estimates, both in rain volume and isohyetal pattern. 

This paper will confine itself to the discussion of daily rainfall, and will examine the funda- 
mental unit of estimation in the GWT, the half-hourly rain calculation. Following a brief review 


2 



of the GWT, data on GWT-derived cloud statistics are presented for one case study day. A dis- 
cussion of some implications of these data leads us to propose (and quantify) several sirnplfying 
assumptions. These are then applied to IR imagery, producing a daily rain map and resultant com- 
parison to the GWT results on that day. 

2. A BRIEF SUMMARY 'OF THE GWT 

In brief, the steps necessary to produce a rain depth (D^) at grid square (i j) at each time 
interval are: 


1) From a sequence of IR imagery, identify and compute the area (A c ) of every entity 
bounded by a 253K isotherm. 

2) Determine successive values of A c for each entity until terminated by a split, merger, 
mingle (simultaneous split and merger) or until lost through evaporation. Entities 
(clouds) resulting from such interactions are considered new entities. 

3) Determine the maximum areal extent (A ) of the entity and use it to normalize each 
observation of A c during the lifetime of the entity. 

4) Use the curves presented in Fig. 1 to determine the fractional area containing precipita- 
tion (A e /A, n ), where A e is the inferred echo (rain) area. Enter the ordinate of Fig. 1 
based on: 

a) The value of A ni (indicates which diagram to use) 

b) The value of A c /A m 

c) The sign of A A c /A t where A t is the interval between satellite images (typically 
30 min). 
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Figure 1 . Cloud area/echo area relationships for infrared data in the GWT. Curves are subjective 
fits to mean data (x’s) and are stratified by maximum cloud area. (Adapted from 
Griffith et ai, 1978). 
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The empirical relationships in Fig, 1 are based on 281 observations of radar echoes and IR 

defined clouds over southern Florida (Griffith et a!., 1978), 

5) Multiply the fractions A £ ./A m by A m to determine, the life history of A e . Determine 
the maximum areal extent of echo (A e ^, nax p and use it to normalize each A e value. 

6) From the value of A e /A e ^ max > and from the sign of A A e /A t enter the ordinate of Fig. 

2 and determine the rainrate (1). This relationship was derived from digitized WSR-57 
radar data from Miami (Woodley etal., 1980, Griffith etal, 1980). The singular point 
at 20 mm If 1 in Fig. 2 is for echoes at their maximum area. 

3 

7) Compute the term 2 ah, for each entity at each time where: 

i=i * 1 

a, is the fraction of A e between 253 and 225K; 
a 2 is the fraction of A c between 224 and 202K; 
a 3 is the fraction of A c less than 201 K; 

and where bj , b 2 , b 3 are the empirically derived coefficients 1 .00, 2.19, and 3.24. This 
temperature weighting term would be 1 .0 (3.24) for a cloud composed entirely of pixels 
at 253K (201K). It is designed to increase the rain volume for clouds with colder internal 
temperature structure. These coefficients, empirically derived, were originally based on 
comparisons between radar echoes and the visible brightness contours in hurricanes 
(Griffith el al, 1976). The infrared coefficients are described in the appendix of Griffith 
etal., (1978). 

8) Compute the rain volume (R y ) for eacli entity at each time: 

Rv“IA e Z ajbjAt 

i“ 1 


( 1 ) 



0.25 0.50 0.75 1.00 0.75 0.50 0.25 

INCREASING DECREASING 

^ e / ^ e (max) 


The echo area/ranrate relationship in the GWT. normalized by maximum echo area. 
The singular point at 20 mm h" ! is for echoes at their maximum area. (Adapted from 
Griffith et ai, 1980). 



9) Apportion half the rain volume to the coldest 10% of the cloud entity area, and the 
remaining half to the next warmest 40% of the cloud. For the GATE data, rainfall was 
apportioned to the whole cloud (Woodley et al, 1 980). For estimation in the U.S. High 
Plains rain was apportioned to the inferred echo area, typically 6% of the cloud area 
(Griffith et al, 1981). For tropical systems, the “10-50/40-50” apportionment is used 
(Augustine et al, 1981). Within these 10% and 40% areas the rain is apportioned based 
on the summation of a parameter b. These summations [(2b) 10 ™ and (2b ) 40 %] run 
over the cloud pixels in the coldest 10% and next warmest 40% cloud areas respectively. 
The b values are given by: 

b » exp(i .784095 - 0.03094 T)/l 1 .1 249 (2a) 

when -3 1 < T < -20° C 

b -- exp(2. 278682 -- 0.01494 T)/ 1 1 .1249 (2b) 

when T < -32° C 

Equations (2a) and (2b) are plotted as the solid line in Fig. 3 indicating a near linear 
relationship between b and T (dashed line). These b values are related to the bj, b 2 , 
b 3 of the preceding section. However, there is a small discrepancy when b 2 and b 3 are 
calculated using Eq. 2b (C. Griffith, personal communication). The b values are ex- 
pressed in terms of percent of gray scale (P) in Griffith et al (1978) and in terms of 
digital count (D) in Griffith et al (1980). It is unclear how many data points were 
actually used to empirically derive these b values. 

10) Apportion a rain depth (Dy) to grid square area (ay) on the basis of grid square tempera- 



0.0338T + 0.222 
legrees C) 



b = exp(1 .784059 - 0.03094D/1 
- 3KK — 20°C 
b = exp(2.278682 - 0.01 494T)/1 

T< - 32 °C 
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Figure 3. Weighting coefficients as a function of temperature in the GWT (solid curve) adapted from Griffith et al„ (1980). 
The dashed line shows a linear approximation. 



where by is the result of applying (2a) or (2b) to Ty and where JJbj, = |>2bJ j o%(40Vf ) 
if by corresponds to a temperature in the coldest 1 0% (next 40% warmest) area of the 
cloud entity. 

3. DATA 

The case study data set consists of GOES infrared digital imagery on 3 1 July 1 980, one of 5 1 
days during FACE-2 on which the GWT was applied (Meitin et ah, 1981). The day was chosen 
because it was the only day of the 51 on which processed, navigated GOES images previously 
existed on the Atmospheric and Oceanographic Information Processing System (AOIPS) at the 
Goddard Space Flight Center. Imagery over the Florida area between 1600 and 0030 (all times are 
GMT) are displayed in Figs. 4-6. The small, trapezoidal area is the FACE target region. The data 
have been enhanced to highlight CTT regions which affect the calculation of rain volume in the 
GWT: black: 252-225K; gray: 224-2.Q2k; and white: <20!K. Every black-outlined feature 
(< 253K) would be defined as a cloud “entity” in the GWT. 

The day appears to be representative of convective development over the Florida peninsula, 
i.e. mostly clear at 1 600 with convective clouds growing during the afternoon. By 0030 the anvils 
of many thunderstorms have merged to form one cloud entity (as defined by the 253K isotherm). 
This sequence of imagery will be frequently referred to in subsequent sections in which we examine 
in detail the impact of the various terms in the GWT equations on the final rainfall estimate. 

4. CLOUD STATISTICS FOR 31 JULY 1980 

It was not possible (nor was it the intent) to replicate the cloud isolation and tracking soft- 
ware utilized by the GWT. Rather, the imagery in Figs. 4-6 are viewed in time sequence. A subjec- 
tive determination is made as to whether a cloud entity (as defined by the 253K isotherm) either 
existed on the prior image, is an entity new to that image, or is the result of a split (or merger) of 
the boundaries of existing entities. We arc attempting to simulate the life history relationships 
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in the GWT so that the dominant parameters in the rain estimation can be identified. Figures 7-13 
summarize, for one day and for one limited area, the cloud-derived statistics for the GWT param- 
eters. 

a) Cloud definition 

Figure 7 shows the number of GWT defined clouds as a function of time of day. Total clouds 
peaked near midday (1800) al ter which the anvils of convective clouds begin to merge, and the 
total number of clouds decreases to 3-4 after 2200. At any given time, only about half the clouds 
are defined as the same entity on the previous image (“old” clouds). Concurrently, very few “new” 
clouds formed; about half the entities at any time are the result of mergers or splits of previously 
existing entities. 

b) Cloud lifetime 

Figure 8 summarizes the observations of S3 entities betv/een 1600 and 0100 GMT. 8 5% had 
lifetimes of 1 h (two images) or less. Only one entity was observed to exist without interaction 
through six images. This storm (in Fig. 4, labelled “A”) will be discussed at length in a later section. 

c) Cloud area 

The 53 entities provided 95 observation-times of cloud parameters. Figure 9 reveals that 
while 79% of the sampled clouds were small (less than 5000 km 2 ) a wide variation existed in cloud 
area. Not all the 95 points can be described as thunderstorms; an examination of the imagery 
shows that the small entities tend to be cirrus debris, white the largest ones tend to be aggregates 
of many thunderstorms in various phases of their life cycles. 

d ) Fractional area of precipitation 

The inferred echo area term (A e ) in (1) is actually a product of two terms, the maximum cloud 
area (A m ) and the fractional area of precipitation (A e /A ). This empirically derived ratio is shown 
in Fig. 10. Because the sampled clouds had brief, GWT -defined lifetimes, entities tended to be near 
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Figure 7. Number of GWT defined clouds as a function of time of day for the data in Figs. 4-6. 
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Figure 9. Histogram of cloud area (< 253K) for the 95 observations of the 53 cloud entities. 
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Figure 1 0. Histogram of the ratio (echo area/maximum cloud area) for the S3 observations. 
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or at their maximum area. As a result, the value of 0.02 for A e /A was invoked most often (see 
Fig. 1 when A n < 2000 km 2 and when A c /A m - 1). This value of 0.02 is used for the calculation 
of R v and should not be confused with the 50% area of the cloud over which the rain is ultimately 
apportioned. The ratio A c /A m had a mean and standard deviation of 0.05. 

e) Rainrate 

Figure 1 1 shows the histogram for the nine possible values of rainrate in the GWT (though 
two are identical). 61% of the data sample had the value (20.7 mm If 1 ) that occurs when A e - 
A e (ma X ) (Fig. 2). The parameter had. a small standard deviation with respect to its mean. 

f) Temperature weighting 

Given two clouds of equal area, the summation term in (1) is designed to increase the rain 
volume for the cloud with the colder internal temperature structure. The distribution of this term 
is shown in Fig. 12. 55% of the sample had no temperatures below 225'K, hence the term was 
identically 1.0. 

g) Rain volume 

The relationship between GWT computed rain volum e and cloud area (Fig. 13) is interesting 
because A c does not explicitly appear in (1); it appears indirectly in A and in (he life history 
relationships from which the A e /A m ratios are derived. Figure 13 is a plot of R v vs. A c for the 95 
observations. The two are highly correlated (r = 0.88) despite the two obvious classes of data. 

The data were stratified by entities increasing in area (27), decreasing in area (3), entities at their 
maximum area as part of a life history (43), and entities for which only a single observation existed 
(22). The GWT empirical relationships for both rainrate and fractional echo area are asymmetric 
with respect to A m (refer to Figs. 1 and 2). Hence estimated R can be one or two orders of 
magnitude greater for clouds defined as increasing in area. Very few entities decreased in area; 
by this stage in their life cycle they tended to split or merge. Some implications of this relationship 
are discussed in (he following section. 



NUMBER (N) 



60 
50 
40 
30 
20 
10 

13.3 17.3 21.1 23.8 20.7 21.1 16.7 11.9 8.2 

INC DEC 

RAIN RATE (mm h -1 ) 

Figure 1 i . Histogram of tne rainrate term for the 95 observations. 
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Figure !3. Log-log plot of GWT computed rain volume versus cloud area. Data are stratified by the tendency of area change. 
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5. DISCUSSION 

Though the data presented in Figs. 7-13 are for one clay, wc believe that, the convective devel- 
opment on this day is representative of the Florida area (where the GWT relationships were devel- 
oped). We believe that the statistics would be similar on any convectiveiy active Florida day. The 
impression one gains, based on these data, is that the life history approach may not be necessary to 
produce daily maps of satellite derived precipitation. This is not to say that clouds (thunderstorms) 
do not have a life cycle similar to that envisioned by Griffith et al. (1976) not is it to say that these 
thunderstorm-rain relationships are independent of storm life cycle. Instead we are stating that 
the conceptual life cycle mode! proposed by the GWT is not applicable when clouds are defined by 
the 253K isotherm. It was noted in this limited data sample that the “average” cloud only lasted 
1 h before interacting with another. Only 3 of the 95 entities ever decreased in area. Most clouds 
life histories were terminated prematurely by mergers or splits. In summary, cloud entities defined 
by the 253K isotherm did not undergo life cycles as depicted by the GWT conceptual model. In 
fact, Woodley et a!. (1980) state that few simple, single clouds are ever encountered. 

The most long-lived cloud entity was the aforementioned storm west of the FACE target area 
between 1600 and 1830 (storm “A” in Fig. 4). Its life history (see Fig. 14) consisted of a steady 
growth over 3 li (dashed curve), terminated by a split. The bracketed numbers above the curve are 

3 

the values of Z ajbj at each time. Normalizing each A c by A m (5862 km 2 ) one can derive the life 
history of echo area (solid lower curve) and from that, infer rainrate (1). The product of I, A e , 

3 

the summation term ( Z a^b,), and A t produce the rain volume (top curve). The numbers in 

is | II 

parentheses above this curve arc the entity’s minimum IR temperature. The absolute minimum 
temperature (T mjn ) achieved by this storm (212K) corresponds to its maximum volume rain 
( 1 2 66 1 km 2 -mn). A similar relationship was found by Negri and Adler { 1 98 1 ) for 1 5 midwest 
thunderstorms. For this isolated storm, it seems that the GWT is able to simulate the evolution of 
thunderstorm rain volume. It is instructive however, to consider what the R curve would have 
been had the cloud been terminated at an earlier stage by a merger. In the worst case, we assumed 
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Figure 14. Sample GWT calculation for storm “A”, west of the target area in Fig. 4. Cloud area 
.. (dashed line) is the primary measurement. 
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that the cloud existed for only one lime step before it merged (at each time step) witli a hypo- 
thetical cloud whose area was insignificant compared to that of the original storm. Applying the 
GWT relationships to these individual entities at six different times generated a radically different 
profile (Fig. 15). The shaded area represents the difference in rain volume solely due to termina- 
tion criteria. While it is clear that growing clouds produce more rain, in the GWT two clouds of 
equal size can have two orders of magnitude difference in rain volume depending on whether they 
happen to grow in clear or cloudy areas. Clouds growing in close proximity to others clouds stand 
a greater chance of anvil interaction (hence redefinition and less estimated rain) than do identical, 
isolated clouds. This result appears to be arbitrary and without physical basis. Note also that this 
cloud does decrease in area (cloud K in Fig. 5). Its life history has been artificially terminated by 
a split at 1900 GMT. 

Another example of a storm whose life cycle is prematurely terminated by the GWT criteria 
is the storm labelled “B” beginning at. 1800 (Fig. 4). Subjectively, the entity can be followed for 
over four hours. It can be seen to reach its maximum extent at 2000 after whMi it decreases in 
area and disappears (warms) at 2230. In the GWT definition, this storm frequently merges with 
the anvils of neighboring storms, becoming a small part of the large cloud mass at 2200. 

Among the cloud parameters in the GWT, the distributions we calculated for rainrate, frac- 
tional echo area, and the temperature weighting term tended to concentrate at one particular value. 
When combined with the data that suggest that cloud entities are not long-lived, (hence dose to 
their maximum area at any time) it is not surprising that R„ is well related (Fig. 13) to A ( . (A = 

A nl for many of the observations). A linear relationship in A . accounts for IT'/t of the variance 
in GWT rain volume, That is, the simple occurrence of cloud, more than any other parameter, 
dominates the calculation of R . Simple relationships between cloudiness • id rain have been found 
ir. several previous studies. Stout et at. (1979) showed correlations as 'nigh as 0.90 between cloud 
area and volume rainrate for sunjectivvly defined clouds. Arktn ( 1979) calculated 6 h averages of 
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Figure 15. Extreme values of GVVT-computeii rain volume for the cloud entity of Fig. 14. Stippled 
regions represent differences in rain volume between the life history computation (top 
curve) and terminating the entity every 30 min (bottom curve). 
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cloud cover collier than 24 5 K and found that to be highly correlated (0.88 ) with 6 h accumulated 
rainfall. On monthly time scales, the simple occurrence of cloud is correlated (0.75 ) with total rain 
(Kilonsky and Ramage, 1976). Similarly, Lovejoy and Austin (1979) found that GOES IR data are 
good for determining rain area but poor for determining rain rates. Secondarily, Fig. 13 suggests 
that knowing the tendence of A c may further, refine the rain volume estimate. Stout et al. (1979) 
also noted the lesser importance (a factor of about two) of the area change term. 

6. PROPOSED SIMPLIFICATIONS 

We consider a cloud with area A c defined by the 253K isotherm and with temperatures T ^ r / o 
and T S()% defining its coldest 1.0% anil next warmest 40% areas iespectively. Assume also that the 
cloud has a minimum temperature T ; . We may compute from (2a) or (2b) the corresponding 
values bjQ^, b 50 ^ , and b max . if we assume that the temperature distribution within tins area is 
uniform (all the Ty and hence all the by are identical) then 

(2b) J0% = -10Nj by (4) 

where N- r is the number of cloudy «253K) pixels. Substituting (4) into (3) and noting that ay 
X M i- (area/pixel X number of pixels) equals A c , (3) becomes: 

D 10% = SRvMc < 5a ) 

This states that half the rain volume is distributed equally among the pixels (grid squares) of the 
coldest 10% area. To test this assumption of temperature homogeneity in the 10% area we looked 
at the “worst” case for each of the 95 entities: at T min the difference (b nnx - U ; o%) is maxi- 
mized, Figure 1 6 shows this term expressed as a percentage of hjQ^. For more than half the 
sample, the point of minimum temperature was the 1 0% temperature, so that b max = b J0 -^. The 
mean value was 2.2%, the highest 18%. These are the extreme values; for temperatures wanner 
than T | the error in the approximation is less. 
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Figure 16 . Histogram of the percent error introduced by the assumption that the temperature distribution is uniform within the 
coldest 10% cloud area. 
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A similar approach may he applied to the next warmest 40 % area: 

D 40 % = 1 ' 251 V A e 

The error in the approximation (2b) 40 ^ = 0.4N T bjj is shown in Fig. 17. A point with temperature 
just warmer than Tjq% has a value b ( j%. The “warmest” value is b^. The error in this approxi- 
mation is greater : mean 13.6%, extreme 44%, However, this is the worst case. As the temperature 
of' the pixel approaches Tj 0 ^ the error in the approximation decreases, 

Figure 18 shows that the distribution of Tjq^, for the 95 entities is quite uniform over a wide 
temperature range. It is here that the GWT’s greatest strength appears to be: by apportioning half 
the rain volume to this area (whether equitably or as a function of pixel temperature) the GWT can 
produce rain in small (warm) clouds while points at the same temperature in clouds with colder 
temperatures elsewhere will not produce rain. Figure 1 9 is the distribution of T-q,^ and also dis- 
plays no pronounced peak. 

We make one final assumption before applying the simplifications to actual data. The most 
simple and direct relation between clouds and rain is one where R y « A , and is substantiated by 
the correlation of Fig. 13 and by correlations found by others. In the simplest case, R y = k X A f . 

We have subjectively let k — 1 mra, the calculated slope of a regression line fit to Fig. 13. Then 
(5a) and (5b) become: 

Dy = 5 mm (TjjCT,^) (6a) 

D y = 1.25 mm (T, 0% < T sj <T S0 *) (6b) 

These equations are independent of cloud life history, independent of cloud area, and independent 
of grid square area. 
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Figure 



7. Histogram of the percent error introduced by the assumption that the temperature distribution is uniform within the 
next warmest 40% cloud area. 
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Figure 18. Temperature (K) that defines the coldest 10% cloud area. 
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7. RESULTS 

Implicit in equations (6a) and (6b) are the variables T| 0 % and T ^q,^, the temperatures that 
define the coldest SO % and 506 of the cloud entity, Because we choose not to define individual 
clouds, we looked at histograms of CTT for the 500 by 500 km region of Figs. 4-6. The variation of 
• the threshold temperatures as a function of time of day is shown in Fig. 20. They undergo a diurnal 
variation of as much as 25K. We then applied these thresholds to the FACE extended area, a 
region 8° of latitude by ! 2° of longitude. In this case, the (i j) grid squares represented an area of 
about 13 km 2 . On the AOIPS, equations (6a) and (6b) were impleuted using a lookup table that 
converts counts (temperature) to rain depth. For example, on the 1600 GMT image, pixels warmer 
than 233K receive no rain, pixels between 232 and 223K each receive 1 .25 mm and pixels colder 
than 223K receive 5 mm. Each satellite IR image, with the lookup table applied, represents a half- 
hourly rain estimate map which can be summed for any time period. 

The result (Fig. 21) is an isohyetal map of satellite derived daily rainfall in GOES coordinates, 
with gray shades representing rain depth (mm). This should be compared to Fig. 22, the isohycts 
derived from the time dependent GWT (reproduced from Meitin ct a!., 1981 ). Because clouds (in 
particular cold clouds) are highly correlated with rainfall, it is not surprising that the schemes pro- 
duce similar results. It should be remembered that neither method represents ground truth, but 
rather estimates based on satellite nephanalysis. The two maps agree well for the large, circular 
rain area off the southeast Florida coast, for the three-celled maxima in the southwest quadrant, 
and for the lighter rain regions off the east coast of Florida. The only notable difference is that the 
GWT 40 mrn maximum due west of the target area is underestimated in the simplified method. 

This was caused by the long lived, isolated thunderstorm discussed previously It is not possible to 
say which display is representative of the true rainfall for one particular storm. A comparison is 
possible for the FACT, target area. Daily (1630-0100 GM T) rain volumes for this area were 2.84, 
2,34, and 2.68 X JO 7 m 3 using the GWT, gages, and the simplified equations, respectively. This 
was a day of relatively small rainfall in the FACE area. This comparison indicates that daily rain 
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Figure 21 . Satellite derived rainfall (mm) for the FACE extended area (large trapezoid) for the period 1600-0130 GMT, the 

result of applying equations (6a) and (6b) to a sequence of IR imagery. The coordinate system is that of the eastern 
GOES, and the gray scale is shewn. 
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Figure 22. GWT derived rainfall for the period 1630-0200 GMT. Contours are every 10 mm beginning at zero. (From Meitin 
et at., 1981). 
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estimates comparable to those of the GVVT can be made from a more straightforward technique. 
Further testing on additional days should confirm this finding. 

8. CONCLUSIONS 

It was not the objective of this study, nor was it possible, to evaluate the GWT by lengthy 
statistical analyses of many data sets. Rather for one day and for one limited area, the magnitude 
and variability of the terms in the GWT equations were explored. The day was typical of convec- 
tive development over Florida, where the GWT was developed. The concept of life history was 
examined via specific examples of satellite imagery (heretofore unpresented in the GWT literature). 

Due to the definition of cloud entities by an anvil edge temperature (253 K), a simple model 
of cloud growth and decay (Griffith et al„ 1 978) did not exist for most of the GWT-defined 
entities. Only about half the entities were defined on two successive images; about half the entities 
at any time were the result of mergers or splits. Because such entities were close to, or at, their 
defined maximum area for most of their lifetime, the empirically derived ratio Ag/A was the 
singular value 0.02 for 43% of the data sample. Similarly the derived quantity A e (echo area) 
typically had link life history, forcing the rainrate term (I) to be 20.7 mm If 1 for 61% of the 
sample. The summation term was identically 1 .0 for 55% of the sampled clouds. 

One parameter with wide variability is the cloud area (A c ). While this parameter does not 
explicitly appear in the GWT equations, it is highly correlated (0.88) with computed rain volume. 
Clouds increasing in area have up to 50 times more GWT rain volume than identical static clouds. 
This effect was more pronounced in small « 5000 km 2 ) clouds. While the concept that growing 
clouds produce more rain is physically realistic, the GVVT has a definite bias towards growing, 
isolated clouds. Such clouds do not interact with other clouds, (lienee being artificially redefined) 
unlike clouds growing in proximity to others. These latter clouds are often defined as having no 
growth because of continual redefinition. 
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Any objective scheme for defining clouds will have shortcomings. However, interactions 
between entities defined at 253K often involve cirrus anvil debris rather than active convection. 
Clouds defined by 253K may often represent entities larger than individual thunderstorms. Such 
entities may be composed of many thunderstorms in various stages of their life cycles. This leads 
us to conclude that the GWT life cycle conceptual model is not generally applicable to clouds 
defined by the 253K isotherm, and that the GWT apportionment scheme is unnecessarily compli- 
cated for daily rain estimation over large ( 1 0 s -1 0 6 km 2 ) areas. On smaller time scales (1-3 h), 
where a life history approach to thunderstorm rainfall would appear to be realistic, the technique is 
severely limited by the 0.5-1 .0 h average lifetime of a GWT-defined entity. 

An important, discriminating parameter in the GWT was found to be the temperature that 
defines the coldest 10% cloud area. (Note that this is not dependent on cloud life history). How- 
ever, it was found that within these areas the temperature structure is fairly uniform. This thresh- 
old (and a similar one for the next warmest 40% area) decreased by 20K as the convection devel- 
oped over the Florida area. A simplified algorithm was derived from the GWT equations by 
assuming: 

1 ) The rain volume is a linear function of cloud area at any time; 

2) The temperature distribution in both the 10% and 40% areas were uniform; 

3) The temperature thresholds varied with time of day and could be applied regionally rather 
than to individual clouds. 

The resultant equations are independent of cloud life history, cloud area, and grid square area. 
They are dependent on the grid square being in one of two temperature regimes. By applying the 
equations to the IR data we derived an isohyetal map very similar to that derived by the GWT A 
comparison of daily rain volumes for the FACE target area yielded estimates of 2.84, 2.34, and 
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